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Abstract 

In the Minimal Supersymmetric extension of the Standard Model (MSSM) squarks from 
the third generation, i.e., the bottom and the top squarks, have a common set of param- 
eters that determine their masses and couplings. These are in the form of the common 
soft mass for the left handed squark from the third generation (m^ ) , the supersymmetry 
conserving Higgsino mass parameter ^ and tan /3, the ratio of vacuum expectation values of 
the two Higgs doublets. This leads to an interesting possibility that the systems involving 
the bottom and the top squarks might get correlated in a non-trivial way even in an uncon- 
strained setup. In this work, phenomenology of the bottom and the top squarks is studied 
at the Large Hadron Collider which exploits such a possibility with a particular emphasis 
on bottom squark decaying to top squark and W boson. Possibility of reconstructing mul- 
tiple W bosons in the final state is identified to be the key in probing the squark mixing 
angles. Further entanglement of the electroweak gaugino-higgsino sector is also not impos- 
sible in scenarios based on phenomenological MSSM while for highly constrained scenarios 
like minimal supergravity/constrained MSSM such entanglements are trivially forced upon. 
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1 Introduction 



The Large Hadron Collider (LHC), by entering into the data-taking phase, has ushered 
a new era in high energy physics phenomenology. As of now, each of ATLAS and CMS 
experiments has collected more than 5.5 fb^^ of data. Data analysed so far, either published 
or preliminary in nature and ranging from 35 pb~^ to about 2 fb^^, are yet to point out any 
clear evidence of physics beyond the Standard Model (SM) of particle physics. Thus, the 
results presented till date are mostly in terms of setting exclusion limits on the otherwise 
viable regions of the parameter space of different plausible scenarios beyond the SM. For 
Supersymmetry (SUSY) as one such viable frameworks to go beyond the SM, the recent 
LHC analyses PElEliaiaElElElElIIOlIIIlIiailSlIIlIiailSlIIS have aggressively 

extended the region already excluded by Tevatron. 

While these are definite and quick improvements over what we knew before the LHC 
data was available, the actual implications of them, at this juncture, should be put in 
perspective. The single most important thing the LHC experiments (the ATLAS and the 
CMS in particular, which are in context of the present work) opened up before us is that 
how powerful and efficient these are to reach out to new regimes in the energy scales. On 
the other hand, it is appreciated for quite sometime now that a particular SUSY scenario 
like the minimal Supergravity (mSUGRA) (or, for that matter, the constrained version of 
the Minimal Supersymmetric Standard Model (the so-called CMSSM)) in which Tevatron 
carried out most of its analyses and LHC might be sticking to for some more time to come, 
is at best an extremely special but nonetheless, a useful benchmark for any first study. 
So, time is just right to move on to explore the chances of how much or what aspects of 
an (essentially) unconstrained SUSY scenarios could, in principle, be put in the context of 
experimental observations. Efforts in this direction are already on \20\ I21j. In fact, in some 
of their recent studies, the ATLAS and CMS collaborations had also presented exclusion 
plots in a framework like the much relaxed phenomenological MSSM [31 [TJ El [9l [16] and also 
referred to the so-called simplified models l2l[5l[TT] following Refs.p2 l [23l [M] . 

Under the circumstances, the SUSY partners of the third generation quarks (the scalar 
quark or squarks) have got a special standing from a few different but somewhat compelling 
considerations. For example, it is well known that the squarks from the third generation, 
viz., the top squark(s) or the 'stop' (t) and the bottom squarks(s) or the sbottom (6), may 
be rather light (compared to their first two generation peers) due to two primary reasons. 
First, in a high scale scenario where soft SUSY breaking masses are defined at a high scale, 
the ones for the third generation run down to a much lower value at the weak scale due 
to the large Yukawa couplings they have. Thus, at the weak scale, the diagonal entries 
in the scalar mass-squared matrices are generically smaller than those from the first two 
generations. This is the so-called inverted mass hierarchy [251 [26l [271 [28] where the scalar 
(SUSY) partners of the massive SM fermions turn out to have smaller soft-masses and vice 
versa. In addition, as electroweak symmetry breaks, large mixing of the chiral states driven 
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primarily by the mass of the partner fermions (the top and the bottom quarks) contributes 
to further lowering of the mass eigenvalues for the third generation sfermions. 

Adhering to a somewhat conservative interpretation of the LHC constraints, thus ex- 
pecting, in general, super- TeV masses for the gluino and the squarks from the first two 
generations |29j , a legitimate search strategy may be to focus on the squarks from the third 
generation |30L I31j. In any case, irrespective of what lower limits are set by an experiment 
on the first two generation squark masses, it is always reasonable to expect them to be 
rather heavy to satisfy different SUSY flavour constraints while allowing for lighter third 
generation squarks with which the constraints are still much relaxed. Also, recent dedicated 
studies at ATLAS and CMS on final states involving 6-jets [3 El EJ HH |T9] in scenarios, 
where only bi and/or ti are lighter than the gluino, tend to allow for a top or bottom 
squark much lighter than 500 GeV while for the gluino and the squarks from the first two 
generations the lower bounds are almost touching 500 GeV and 1 TeV, respectively. 

Thus, at a time when dedicated search for the squarks from the third generation would 
be on, while other strongly interacting superpartners still remaining unreachable, the pre- 
paredness to decipher the most from a positive signal would be something of genuine priority. 
There have been some studies in the past |32^ [33] where techniques of measuring masses 
of the squarks from the third generation were discussed in scenarios where these squarks 
are produced in the decays of gluino. In particular, these studies suggested studying the 
edge structure of the nitb distribution when gluino decays to tbxi^ through a light stop or 
sbottom. Further, in recent times, phenomenology of stop as the next to lightest supersym- 
metric particle (NLSP) has been considered at the LHC and/or Tevatron in scenarios like 
Gauge Mediated SUSY Breaking (GMSB) iMj or the CMSSM [35]. Interesting studies of 
probing the third generation squark sector, in particular, the stop sector through the study 
of the Higgs bosons were taken up in Refs. [36l EZj- 

Once we know about the primary quantities like the production cross sections of these 
SUSY particles and their masses, we would turn to learn about their couplings. These, 
together, would be reflective of the amount of chiral admixtures present in the mass eigen- 
states of the bottom and the top squarks. Dedicated study in this area is not abound, 
particularly in the context of LHC. However, a rather recent work [38] put the issue in 
perspective with the top squark sector in reference. Further, there are recent works \39\ HO] 
proposing new techniques to reconstruct the top squark by tagging the top quarks. 

As we will see in the next section, the observables in the sbottom and the stop sectors, 
particularly the basic ones like the masses and the mixings, are controlled by some param- 
eters which are common to both. Incidentally, some of them also control the compositions 
of the charginos and the neutralinos in an important way. Thus, it is only natural to expect 
some reasonable correlations inherent in these sectors. Such correlations might provide 
some phenomenological handle in systematic explorations of these sectors in tandem and 
this is what we like to address here. 

In this work, we perform a case-study with the production of the lighter bottom squark 
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at the 14 TeV LHC which eventually decays either to a top squark and a SM W boson or to 
a lighter chargino and a top quark or to a bottom quark along with a neutralino. Note that 
the first of these decay modes is the only one that involves squarks from both sbottom and 
the stop sectors and thus could reflect on the simultaneous features of these two sectors, i.e., 
the so-called entanglement. Note that the decay bi — )■ tiW^ boson would be significantly 
enhanced only if both of the squarks have significant SU{2) (left) admixtures. This decay 
mode was first discussed in Ref. |41j. In this work, we point out that this is where the 
entanglement between the sbottom and the stop sectors is likely to play an important role. 

For all the decay modes mentioned above, the resulting final states would contain 6-jets; 
either from the direct decays of the bottom squark or from the subsequent decays of the 
top squark or the top quark. It was also noted in Ref. [41J that the pT spectrum of the b 
quarks, the VF^-bosons and the missing px spectra could be different for the above decay 
modes. Thus, the count of events enriched with bottom quarks, bosons and missing px 
with characteristic kinematic properties could be indicative of the mixing-pattern present 
in the sbottom and stop sectors. We also point out that the LHC running at 7 TeV is not 
likely to offer much insight in such a study. 

he paper is organised as follows. In section 2, we briefly outline the roles of different 
SUSY parameters that control the masses and the mixings in the sbottom and the stop 
sectors. We indicate how some of these parameters affect the compositions of the charginos 
and the neutralinos. We also touch upon the kind of correlation these common parameters 
might bring into the combined system. In section 3, we discuss the signals and the potential 
backgrounds at the LHC. A brief account of the setup and our analysis followed by the 
important observations are presented in section 4. In section 5, we summarise. 



2 Masses and mixings 

In this section, we outline the correlated nature of the system comprised of the third gen- 
eration squarks and, to a varied extent, the charginos and the neutralinos. 

The mass-squared matrix for the squarks, more relevantly, those from the third gener- 
ation has the following generic form in the (m^ , rriq ) basis: 

o / nij T + mi ma mrn \ , , 

K={ ' 2^ 2 

where 

"^LL = -|- (T3 — e sin^ 6*^^) m| cos 2/3 

m\^ = w-lg^ + e sin^ 9w "t-z cos 2/3 

^Ir = Ag^-fJ'R (2) 

where m^^ is the bottom or the top quark mass, m^ (m-g^^) is the SUSY-breaking soft 
mass for the left (right)-chiral bottom or the top squark, T3 and e are the isospin quantum 
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number and the electric charge of the squark in context, 9w is the Weinberg angle and Aq^ 
is the soft SUSY breaking scalar trilinear parameter that appears in the scalar potential of 
the MSSM {Ah or At, as the case may be). R is equal to tan/3 (cot /3) for the sbottom (stop) 
sector. The second term in each of the first two expressions of equation (2) is proportional 
to cos 2/3 and are known as the SU{2)l and U{1)y .D-terms originating in the quartic 
term of the scalar potential. 

Further, let us write down the general formula for the mass eigenvalues of the physical 
states obtained on mixing of the chiral degrees of freedom: 



mir r. =\ {2mlt + uiIl + tuIr T ^ (ml ^ - m|^)2 + Aml^ml^ (3) 



where the subscript 1(2) on the left hand side corresponds to the negative (positive) sign 
in front of the term under square-root on the right thus ensuring the mass-variable with 
subscript 1 (2) is the lighter (heavier) of the two mass eigenstates. The expressions involving 
the mixing angles are as follows: 



2 2 



bi,ti 



(4) 



The lighter mass eigenstates in terms of the mixing angles have the general form: 



ii = cos L + sin Oiiji 

bi = cos 0~J)L + sin e-J)R (5) 

while the heavier ones are orthogonal states, respectively. Note that < 0^ ^ < tt by 
convention. Also, under this convention, a mass eigenstate is purely left-chiral in nature 
for = 0° and entirely right-chiral when 0^^ = 90°. A detailed description of this sector 
can be found in the thesis work of Ref. [l2] . 

A key issue which we like to exploit here in this work is that the soft mass for the 
left-chiral degrees of freedom (mp, ) is common to both sbottom and stop sectors. The 
sole difference between the two sectors in this part arises due to the so-called D-term 
contributions which are different for the two members of the isospin-doublet {m^^ and m^^ ) . 
On the other hand, n-g^^ for the two sectors can be different {"m^^ and "mi^, respectively) 
as these correspond to isospin singlets not related by any symmetry. 

A few important points pertinent to the present study that emerge from equations 
(1), (2) and (3) are as follows. Clearly, through its presence in the off-diagonal terms of 
the mass-squared matrices, /i controls the mixings, i.e., the chiral contents of the mass 
eigenstates, provided the trilinear couplings are not too large. Also, note that fi appears 
as a product with tan/3 (in the sbottom sector) or cot/3 (in the stop sector). Hence, choice 
of tan /3 as a basic parameter could determine the role of fi by tempering the product-term 
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up to an order of magnitude over the range 5 < tan /3 < 50. The patterns of such mixings 
are effectively captured, in a somewhat global way, in the parametrisation of Ref.|43j in 
terms of a ratio of the diagonal and the off-diagonal entries in the mass-square matrix. 
It can be easily checked that for At (Ah) « ^cot /3 (^tan/3), the off-diagonal terms in 
the two sectors are comparable when tan /3 ~ 6 which is uniquely determined by the ratio 
mt/mb- Thus, in a situation where the diagonal terms of the mass matrices for these two 
sectors are comparable, the role of n in determining the mixing would be similar under such 
circumstances, for any value of /x. 

Further, it is to be noted that the relative values of the soft masses and nin 

appearing in the diagonal terms of the mass-squared matrix influence the chiral contents of 
the sbottom and the stop mass eigenstates in a definitive way in the presence of mixing. For 
example, maximal mixing in the top (bottom) squark sector requires m^^ibL) ~ m^^(m^^) 
at the weak scale. This is a robust, but rather intuitive, requirement. Larger the off-diagonal 
term is, larger a deviation from the equality of these masses can be afforded for obtaining 
maximal mixing. Thus, assuming that the off-diagonal terms of the mass-squared matrices 
in both the sectors are likely to be primarily controlled by the corresponding fermion masses, 
almost near equality of the soft chiral masses in the sbottom sector is required to yield close 
to maximal mixing. By the same token, for the stop sector, a somewhat larger difference in 
the chiral soft masses can still result in a large mixing for given values of A^^t, /i and tan /3. 

However, such near equalities are unlikely to be achieved in generic high scale scenarios 
with universal scalar masses (including mSUGRA) and with a grand desert between the 
high scale and the electroweak scale. Take for example, an S'O(IO) GUT inspired scenario 
with a universal scalar mass mo |44j where the 50(10) breaks directly to the SM gauge 
group at the GUT scale. Both left and right top squarks leave in the same 5-plet of the 
SU{5) embedded in 50(10). This would lead to (up to effects from physics above the GUT 
scale) = m^^ at the GUT scale. It immediately follows that < at the weak 
scale due to renormalisation group (RG) running. In contrast, m^^^ can be different from 
m^^ at the high scale since bi and reside in different multiplets of 50(10). Hence, at 
the weak scale, m^^ and m^^^ can have any relative hierarchy if nonuniversality in scalar 
masses at some high scale is allowed for. Thus, in particular GUT motivated scenarios, the 
chiral contents of the sbottom mass eigenstates may vary in a more relaxed fashion unlike 
those of the stop mass eigenstates. 

Thus, while maximal mixing in the sbottom sector can somewhat more naturally be 
achieved by fulfilling the requirement m^^ ~ ' same in the top squark sector can 
hardly be realized in popular GUT motivated scenarios that assume universal scalar masses 
at the high (unification) scale and the presence of a 'grand desert' between the latter and 
the weak scale. In other words, observation of maximal mixing in the stop sector could 
very well indicate a departure from this popular paradigm. In this work, we try to exploit 
this issue. Thus, we adopt a purely phenomenological approach. We vary the chiral soft 
masses in the stop sector freely to explore the mixing in this sector and its implications. 
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On the other side of the story, the values of the soft masses of the U (1) and SU (2) gaug- 
inos (Ml and M2), with respect to fi, are crucial for the masses and the compositions of the 
charginos and the neutralinos (their gaugino and higgsino contents) which, in turn, broadly 
dictate the phenomenology involving them. Note that the charginos and the neutralinos 
interact with sbottom and stop eigenstates with both gaugino and higgsino components 
they have. This immediately hints towards a possible bridge that /i could provide between 
the two sectors. In addition, the role of /i can be seen in conjunction with that of tan /3 in 
the stop and sbottom sector (as explained above) while the role of the latter in the chargino 
and the neutralino sectors could assume importance under specific situations. An illustra- 
tive study to this end encompassing both the sectors is beyond the scope of the present 
work and will be taken up in a subsequent study [55] . 

In this work we restrict ourselves to the study of possible correlations that might be 
present in the sbottom and the stop sectors. Such correlations are likely to be best mani- 
fested in studies where both sectors have explicit involvements. The most significant of such 
situations could be realised in the productions of the lighter sbottom (lighter stop) at the 
colliders followed by its decays to the stop (sbottom) where Pi^^-bosons appear. We thus 
focus on a particular case of production of bottom squarks at the LHC followed by their 
subsequent decays with an aim to study the imprints of mixings in the involved sectors. 



In this section we present some aspects of phenomenology of the lighter of the bottom 
squarks at the LHC. Recently, in the context of a SUSY/MSSM 'golden region' gSl liTj . 
such a study was taken up in Ref. |48j . This was necessarily limited, though in a rather 
motivated way, to a low /i regime with all the squarks from the first two generations taken 
to be rather heavy. The study stuck to an appropriate final state. In the present study, 
for our purpose, we adopt a rather open approach and analyse other possible final states 
originating from sbottom production at the LHC. This, as we will see, would inevitably 
involve the top squarks and thus may potentially shed light into the chiral compositions of 
the bottom and the top squarks in a correlated way. As discussed in the Introduction, this 
might even have reference to the texture of the gaugino sector in the same go. 

Unlike in Ref. |48j where only pair-production of bottom squarks was considered, we 
keep the options open for other production processes which may eventually lead to a pair of 
bottom squarks at the LHC. These include production of a lighter bottom squark in associ- 
ation with a gluino and pair-production of gluinos where a gluino may decay subsequently 
to a bottom quark and an sbottom. These three modes of lighter sbottom production are 
shown schematically below: 



3 The Signal and the Background 



• pp^ bibl , bibi , blbl 
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• pp-f gg : g ^hh\, hhi 

Now, an sbottom would always lead to a bottom quark in its decay, be it from a direct 
decay or through top quark production under a cascade. These are illustrated below. 

Note that any of these modes could lead to at least a pair of lighter bottom squarks and 
up to four bottom quarks (in the case of gluino pair-production) in the final state. 

As indicated above, the bottom squark may have different possible decay modes out of 
which the ones to iiW^, txf, &Xi ^-nd would be of importance. As for ii, its decays to 
txi, tx2: ^Xi would be in context. Further, the gauge bosons (like W^^), the charginos and 
the neutralinos all would contribute to both leptonic and hadronic final states. In fact, as 
pointed out in the Introduction, the decay channels bi — )■ iiW^ and bi — txf could lead 
to identical final states. Hence, if the branchings in these two channels are complementary 
in nature, rates for the final state events may not be much sensitive to the variation of the 
same. Thus, they may not shed much light on the couplings (and hence, on the mixings 
angles) involved. We will get back to this issue in the next section where we study situations 
under which clearer imprints of branchings and hence mixing angles are left in the events 
rates of the different final states. We would also contrast those to a generic scenario where 
all decay modes of 61 are open and some of them are significant. 

It is clear from the above discussion that the multiplicities of bottom quarks in the final 
state, at the parton-level, may vary between 2 and 4. While an experimental study of final 
states even with a moderate bottom quark multiplicity is a challenging proposition, their 
presence, nonetheless, would boost the counts for low-multiplicity (up to 2) final states 
through enhanced combinatoric factors. With these general possibilities in mind, we pick 
a benchmark (reference) MSSM spectrum in the next section for our subsequent analysis. 

The final states we consider are the following: 

1. 2 b-jets + > 4 jets + > 1-lepton + 

2. 2 b-jets + > 4 jets + same-sign dilepton (SSDL) pair + 

3. 2 b-jets + > 4 jets + opposite-sign dilepton (OSDL) pair + 

where by leptons only the electrons and the muons are meant. The first of these final state 
has been studied in Ref. |48| as mentioned in the beginning of this section. In our case, this 
remains to be an important channel for two reasons: first, because the presence of a lepton 
in the final state helps negotiate the otherwise large pure QCD background and second, 
because it is less suppressed compared to dilepton final states indicated above. However, as 
expected, for the dilepton final states the backgrounds are further suppressed thus making 
them worthy of a closer study. In addition, there is always the well-known advantage of 
studying signals in multiple channels which, when studied simultaneously, may potentially 
shed light on intricate issues pertaining to the mass-spectrum and the involved couplings 
in a more efficient and definitive way. 
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As for the backgrounds for the final states indicated above, we considered several SM 
processes that could be potentially strong. These are Zti+0, 1, 2, 3 jets, Wtt+0, 1, 2, 3 jets, 
tt + 3,4 jets, tbW + 0,1 jets. Note that tt-pair production and the same with extra hard 
jets in low multiplicity (up to 2 extra jets) are unlikely to yield serious background because 
of the minimum lepton and jet multiplicities we required for our signal. Fakes from charm 
and light quark jets are not taken into account in this study. 

4 The Setup and the Analysis 

The setup for our analysis is based on two 'benchmark' scenarios about which variations 
are studied. These consist of two sets of weak-scale MSSM input parameters arranged in 
a minimalistic way for the purpose. These would be sufficient to demonstrate the goals 
of the present study. The first one is shown in Table [l] This explores the prospect of a 
heavier spectrum where the masses of the gluino and that of the squarks from the first two 
generations are somewhere near or just exceeding their present bounds in conformity with 
recent studies at the LHC and their interpretations within the much relaxed 19-dimensional 
pMSSM framework |29j . The second benchmark scenario is elaborated in Table [2| This, in 
turn, represents a lighter spectrum with somewhat lighter squarks from the third generation 
along with lighter gluino, charginos and neutralinos that are still very much allowed by the 
LHC data, as explained later in this section. 

In addition, for these two spectra, we ensure compatibility with other experimental 
constraints like the lower bound on the mass of the lighter chargino obtained from the LEP 
experiments (~ 105 GeV) and those pertaining to the anomalous muon magnetic moment 
(a^ = ^) dg [50] and the rare Flavour Changing Neutral Current (FCNC) process like 
6 — )• S7 {51^ [52] while choosing these benchmark points. In fact, the observations in the 
latter two experiments, when considered in conjunction, favour ^ > [5^. This is why we 
choose to work with positive values of fi in the present study. As for the SM-like lightest 
SUSY Higgs boson, we require a somewhat relaxed lower bound of 111 GeV as against the 
actual LEP constraint of 114.4 GeV [54^ 155] . This relaxed bound is consistent with the 
more precise estimation [56^ 157] of theoretical uncertainty involved in predicting the mass 
of the Higgs bosorj^ 

As indicated in the Introduction, our goal is to work with a somewhat light bottom 
squark which is expected to be within the reach of LHC running at its design centre of 
mass energy, i.e., 14 TeV. Moreover, for our purposes, such a sbottom should have enough 
phase space to decay into a top squark along with a W boson. The probability of such a 

■^Strict compliance to the LEP Higgs bound and other precision observables like and BR[6 — > 57] are 
enforced only to these reference points. Everywhere else, where the purpose is to explore the sensitivity of the 
event rates to the mixing angles, we just generated the combinations of these angles by varying the relevant 
MSSM inputs without worrying about their compatibility with these experimental observations. 
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decay would then be maximised when both 61 and ii have substantial left chiral {SU{2)) 
admixture. However, by requiring such an admixture in them, it is impossible to achieve 
a mass-splitting of Am^^^^ ~ Am^^^^ > mw between these two states. This is because, in 
such a limit, their masses at the weak scale are related, and at the lowest order, differ only 
by an SU{2) D-term which is anyway not large (see equation (2), first expression). Hence, 
a compromise is needed and chiral mixings are to be allowed for. This would then definitely 
suppress the coupling biiiW^ but at the same time open up the required phase space for 
the decay bi — )• iiW^ to take place. 



Parameters 


Input values 


Output spectrum 


Gaugino masses 


Ml = 200 


m° ^ 198 


(in GeV) 


M2 = 400 


my.0 ^ 403 m + ^ 403 




M3 = 1200 


mg ^ 1193 


Chiral Squark masses 






(1st and 2nd generations) 


m-1,2 = 1000 

1l,r 


m-1.2 ^ 1005 - 1010 


(in GeV) 






Third generation 


= 700 


m^^ ^ 709 ^ 1014 


chiral squark masses 


mr = 1000 

OR 


Oi = 1.9° 


(in GeV) 


mr = 500 


m^^ ^ 455 m^^ ~ 757 


and Mixing Angles 




O't = 67° 


Slepton Masses 


m, = 600 




(in GeV) 




m^ ^ 600 


74-parameters 


^ = 




(in GeV) 


At = -800 




niA (in GeV) 


500 


^ 699 


/i (in GeV) 


700 


m^o 711 


tan/3 


10 


m + ^711 

X2 



Table 1: A somewhat heavy benchmark SUSY spectrum and the weak-scale values of the MSSM 
input parameters used to obtain the same. The input soft mass of the CP-odd Higgs is taken to 
be 500 GeV. The SUSY spectrum generator used for the purpose is Suspect v2.31. Throughout 
the analysis m^p = 172.5 GeV is used. Note that m^ < mg. 

In Table [T] we collect the set of relevant MSSM input parameters for the high-mass 
reference point and the resulting output spectrum that conform with the setup described 
above. We call this the high-mass benchmark point/spectrum. Note that the mixing angle 
in the sbottom sector, is rather small thus making the lighter sbottom an almost pure left 
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chiral state. However, 9^ is around 67° and hence it has some (15%) left chiral admixture. 
By varying the relevant MSSM inputs about their benchmark values, the mixing in the top 
squark sector can be significantly altered. This would affect the decay branching fractions 
of the bottom squark. As we would see later in this section, this can have some impact on 
the signatures at the LHC. 



Parameters 


Input values 


Output spectrum 


Gaugino masses 


Ml = 100 


m° ^ 100 


(in GeV) 


M2 = 200 


m^o ^ 207 m ± ^ 207 




M3 = 600 


TJlg ~ 655 


Chiral Squark masses 






(1st and 2nd generations) 


m-1,2 = 1000 

Il.r 


m-1,2 ^ 1005 - 1010 

Il.r 


fin CrpV) 

I 111 VJiC V J 








= 480 


TDr ^ .500 TDr ^ 71 .5 


chiral squark masses 


= 700 


Oi = 4.3° 


(in GeV) 


rrif = 390 


mi ^ 349 mi 571 


and Mixing Angles 




= 57.3° 


Slepton Masses 


m-n = 600 




(in GeV) 


m^- ^ 600 


y4-parameters 


Ab = 




(in GeV) 


At = -500 




niA (in GeV) 


500 


rriyo ^ 845 


/X (in GeV) 


850 


m^o ^ 849 


tan (3 


10 


m^± 850 

X2 



Table 2: Same as in Table [T] except for a somewhat smaller masses for the lighter charginos and 
neutralinos, the gluino and the third generation squarks. Note that niq > nig. 

To decide on the benchmark spectrum with lighter masses (Table [2]) , a closer look at 
recent LHC studies [3 El |9l [19] is warranted. These analyses discuss signals with heavy 
flavour jets (6-jets) and large missing transverse energy without [8] and with |i9jj an isolated 
lepton (e and/or n). There, either bi [S] or [H] are assumed to be the lightest squark 
and the branching fractions of gluino decaying into them (in the respective cases) is 100%. 
Thus, in these studies, 6i-s or ti-s are produced either directly in pairs or via production 
and subsequent decays of gluino. Further, it is assumed in these experimental studies that 
bi and ti always decay in specific channels, i.e., BR[6i — )• bxi\ [8j and BR[ti — )• ^Xi"] IS] are 
100% in the respective cases. Ref.jH] rules out rrig < 720 GeV for < 600 GeV. Ref.[9| 
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excludes rrig < 500 — 520 GeV for m^-^ between 125 GeV and 300 GeV. In Ref.|9j, exclusion 
limit is also placed on the rrig — m^o plane when ti is still the lightest of the squarks but is 
heavy enough (rrig < m^^ + mt) such that g cannot decay into an on-shell ii and it decays 
via ofF-shell ii to ttxi final state with a decay branching fraction of 100%. The analysis 
excluded m^o < 40(80) GeV for rrig < 570(540) GeV. All these reported exclusions are at 
95% confidence level. 

However, the above limits may be considered conservative because they assume the 
Br[^ — )• bib + h.c] or Br[^ — )• tit + h.c] to be 100%. In a situation where there is more than 
one third-generation squark lighter than the gluino, this assumption would not hold and 
may result in relaxed bounds. In addition, note that events from direct stop pair production 
are reported [9] to have a very low acceptance; presumably due to hard kinematic cuts used 
in the analysis. These issues may very well dilute the sensitivities of these experimental 
analyses (to a given data-set) which would eventually lower the exclusion limits. In Table 
[2} we tried to exploit this caveat to our advantage and settled on a somewhat lower mass 
for the gluino (~ 655 GeV) in relation to mj^^ (w 500 GeV) considered. 

For the two benchmark scenarios, we work in an otherwise unconstrained SUSY scenario 
except retaining an imprint of unification of gaugino masses (at a high scale) in the choice 
of their weak scale values. Note, however, that such a choice is only to keep the analysis 
tractable and does not feature an essential part of our study. As is well-known, departure 
from such an assumption, can easily have nontrivial phenomenological implications. Two 
such examples are discussed in the context of the benchmark scenarios, the emphasis being 
on the subtle handles these may provide in the analysis and interpretation of the collider 
signals. A systematic study of the implications of such departures is beyond the scope of 
the present study though and would be taken up elsewhere [15] . 

In Table [3] we collect the lowest order cross sections for the strong-production processes 
that lead to a pair of bi, viz., pp — > bib^, bibi + h.c, big + h.c, gg for the 14 TeV run of the 
LHC. The cross sections are calculated by the event generator Pythia v6.420 [58] and cross- 
checked with CalcHEP v2.5.4 [59]. Note that this is a conservative estimate since the next 
to leading order (NLO) contribution from QCD to squark (including stop and sbottom) 
and gluino productions [60^ l6Tl [62] and the same combined with next to leading log (NLL) 
resummed soft-gluon contribution [63[ l6l] may increase the cross sections by~ 30 — 40% on 
an average, for our benchmark scenarios, the NLL contributions being far more important 
for the case of gluino and squarks from the first two generations. Two sets of cross sections 
are presented for the two benchmark spectra of Table [T] and Table [2} As expected, the 
respective cross sections are much larger for the lighter spectrum of Table [2] Also, in 
both cases the relative magnitudes of g-pair cross section are significant. So, if there is a 
reasonable branching fraction for the decay g — )■ bib + h.c, pair-production of gluino could 
become a useful source of 6i-pair [32l [33] . 

The high-mass scenario (of Table jT| with rriq^ ^ < rrig turns out to be a rather conserva- 
tive one. First and foremost, the individual cross sections for all the contributing processes 
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Spectrum 


(pb) 


(pb) 


^99 

(pb) 


Table [l] 

{nig^ ^ < mg) 


0.030 


0.003 


0.022 


Table |2| 

(m,^^ > m-g) 


0.265 


0.110 


1.764 



Table 3: Lowest order production cross sections at the LHC {^/s = 14 TeV) for different strong 
production processes that lead to pair of lighter bottom squarks using the mass-spectra of Tables 
1 and 2. CTEQ6L parton distribution function is used with renormalisation/factorisation scale 
set at a/S. 

are on the smaller side because the sparticles are heavier. Second, the branching ratio for 
g — >• bib is affected since g — t- Q-^ {+ h.c.) is kinematically accessible and when summed 
over 4 ffavours, would become dominant. Third, we miss out on potential contributions 
from the production q_^^, in pair or in association with a gluino, which could have led up 
to 4 bottom quarks had the hierarchy been < < rrig < niq^ ^ ~ m^^ ~ m^^ . 

The suppression of Br[g — )• bib + h.c] is demonstrated in the upper part of Table [4] for 
the heavy spectrum of Table [TJ It is, however, a straight-forward exercise to find out to 
what extent Br[^ — t- bib + h.c] may get enhanced if the mass-hierarchy of the gluino and 
the squarks from the first two generations is flipped, other parameter remaining the same. 
Note that the ratio of the branching fractions of the gluino to bi and ti (i.e., 9.8 : 13.1) 
reflects the the ratio of the decay widths of the gluino in these two modes. Thus, when only 
these two modes add up to 100% of the branching fraction, Br[^ — )• bib + h.c] would be 
around 43.5%. More importantly, such tweaking of the squark-masses would immediately 
bring in further contributions from heavier squarks when they decay into gluino followed 
by the latter decaying into bi. In contrast, our choice of low-mass spectrum of Table [2] 
represents a favourable scenario on all respective counts laid down above. It is evident from 
the lower parts of Tables [3] and |4] that the low- mass benchmark spectrum boosts both cross 
sections and its branching fractions relevant for our purpose. Thus, the two benchmark 
scenarios are arranged to demonstrate the extremal situations. 

Table [4] also collects the branching fractions of 61 and ii to different modes that are 
instrumental for the study we undertake. For the high-mass spectrum, we see that the 
crucial mode 61 — )• W~ti has a reasonable branching fraction of about 64%. For the low- 
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Spectrum 


Branching fractions 
q and/or ^-decay 
(in %) 


Branching fractions 
6i-decay 
(in %) 


Branching fractions 
ti-decay 
(in %) 


Table m 
(heavier masses) 
mg > niq^^^ > > nii^ 


g -^bbl + III ■■ 9.8 
g til + til ■ 13.1 
g —J- qq* + qq : 77. 1 


bi W'ti 
bi -> txi 
li ^ bxl 
bi bxl 


64.1 
19.6 
15.1 
: 1.2 


k txi ■ 87.0 
ii bxt ■ 13-0 


Table |2| 
(lighter masses) 
rriq^ ^ >mg> m^^ > m^^ 


qR-^ qg : 95.0 
qi^qg: 70.0 
g^bbl+ III ■■ 48.5 
g^ti\ + til ■■ 51.5 


bi W'ii 
bi txi 
li ^ bxl 
li ^ bxl 


34.5 
36.2 
27.7 
: 1.6 


ii txi ■ 39.2 
ii bxt ■ 60.8 



Table 4: Branching fractions of g, bi and ti for the mass-spectra given in Tables 1 and 2. 

mass spectrum, this branching fraction drops to around 34%, presumably due to a combined 
effect of diminished sphtting between bi and ti and possible enhancements of the competing 
modes. This is a plausible explanation since, as can be seen from Tables [T] and [2j the left 
{SU{2)) admixtures in hi and ti are not different in these two cases. 

One of the goals of this work is to understand how sensitive the event rates in differ- 
ent leptonic final states are to the mixing angles in the third generation squark sector. 
As pointed out earlier, since all the relevant leptonic final states get contributions from 
almost all possible branching modes in this sector, our ability to filter out the individual 
contributions to the extent possible would help extract information on the mixing angles 
involved. 

Towards this end one needs to know how sensitive the individual branching fractions of 
the bottom and the top squarks are to the variations of the mixing angles 9^^ and 9^. We 
choose to study the effect of variation of 9f; only. The reason behind this is two-fold. First, 
% is more sensitive to small variations in the SUSY inputs than thanks to the amplifying 
SM factor mt sitting in the off-diagonal term of the corresponding mass-squared matrix. 
Second, 9^ variation directly affects the branching fraction of ti. Nonetheless, an analogous 
study with respect to 9^ would have its own characteristics though effecting a comparable 
variation on it would require major tweaking of the soft parameters of the bottom squark 
sector. 

At this stage, for such a study to be meaningful, it is to be necessarily assumed that the 
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Figure 1: (a) Variations of branching fractions of the hghter bottom squark (red/dark grey band) 
and hghter top squark into different major decay modes as a function of 6^ while keeping m^^ within 
the range 445-465 GeV. The variation of 9^ is achieved (see Figure [2]) by varying the soft parameters 
TTij^ and At while keeping all other parameters relevant to the third generation squarks fixed at the 
benchmark values. 



masses of the lighter sbottom and the lighter stop are already known from the experiments. 
This ensures that the concerned kinematics remain mostly unaltered for the relevant pro- 
duction and decay processes. Thus, we vary m^^ and At to achieve the variation in 6^ while 
allowing m^^ to vary only within a certain range ((±10 GeV) about the reference value of 
m-j^ = 455 GeV (see Table [T|. Clearly, such variations do not at all alter either the mass 
of the lighter bottom squark or the mixing angle in that sector. We also keep fi, tan /3, Mi 
and M2 fixed during this variation so that the masses and the mixings in the chargino and 
the neutralino sectors remain almost fixed in the process. Altogether, this ensures that the 
variations we see in the branching fractions of the bottom and the top squarks are almost 
entirely due to the variation in 9^. This would definitely help obtain a clearer picture of the 
role of 9^ in shaping the pattern of cascade decays of the lighter sbottom. 

In Figure[l| we present the variation of branching fractions of both 61 and ii as functions 
of Note that the products of different branching fractions of these excitations determine 
the predominance of certain cascade-patterns when they decay. In the present case, the 
possible decay mode of lighter sbottom to lighter stop and charged Higgs boson is closed, 
the latter being rather heavy (resulting from our choice of an input mA of 500 GeV). 

Out of the decay widths that enter the calculation of the branching fractions of the lighter 
sbottom, only r(6i — )• iiW~) depends upon 9f: and is actually proportional to cos^ 6^ cos^ 9f 
|41j . Thus, as a function of 9^, the variations of different sbottom branching fractions are 
solely determined by the branching profile of 61 — )• iiW~ which goes as cos'^ 9t (and hence 
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the symmetry about 9^ = 90°), 9^^ being held fixed, as is the case her^ Note that such a 
symmetry with respect to 9^ is not there in the variation of branching fractions of ti. This 
is because the they have a somewhat comphcated dependence on 9^. Clearly, this variation 
neither alters nij^^ nor Thus, bi remains to be almost purely left chiral. 

It is to be noted that the variations in Figure [T] are in the form of bands. The reason 
behind this is that the variation of 9^ is achieved by simultaneous variations of m^^ and At 
over the ranges 300 GeV < m^^ < 1500 GeV and —3 TeV < At < +3 TeV, respectively 
while ensuring m^^ to be roughly in the range 445 < < 465 GeV. 




(in GeV) 

Figure 2: Ranges of 9^ in degrees presented in colour-code (for the setup in Figure [T| as and At 
vary keeping 445 GeV < < 465 GeV. m^^ and are kept fixed at 700 GeV and 1000 GeV, 
respectively. 

Figure [2] reflects some important features relevant for the phenomenology under the 
given setup. These are: 

• The graph is almost symmetric about At = 0. This is expected when At dominates 
in the off-diagonal term of the stop mass-squared matrix. Given that tan /3 = 10 
and n = 700 GeV for our case and the off-diagonal term dominates over most of the 
At-range studied above. 

• The above issue is equivalent to having complementary angles and this is clearly seen 
in the figure. 

''in fact, while the left branch with 6j: < 90° arises for At < 0, the right one with 9j: > 90° results from At > 0. 
Also, moving away on both sides from 6f: = 90°, i.e., increasing amount of mixing, corresponds to increasing 
\At\. Since, with increasing \At\, mixing is dominantly determined by At, it is expected that further one is from 
Of — 90° on either side of the curve, it is more likely that similar values of \At \ result in similar values of branching 
fractions. 
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• The resulting band is rather narrow. This is because we require m^^ be within the 
range 445 GeV < m^^ < 465 GeV. 

• For values of m^^ much larger than w-q^ one requires large values of \At\ for bringing 
down m^^ to the mentioned range. The resulting mixing can be close to maximal or 
lower but never reaches very small values (corresponding to dominance of left-chiral 
admixture) . 

• Maximal mixings (around 45° and 135°) are expected for m^^ in the vicinity of "^g^ • 
However, as mentioned earlier, since the mixing in the stop sector is driven by mt, At 
plays a crucial role in mixing and help achieve close to maximal mixing even when 
m^^ is somewhat away from m-g^ . This is clear from the bluish-purple and yellowish- 
orange bands extended over the range 700 GeV < m^^ <850 GeV. 

• For smaller values of m^^ (> 450 GeV), m^^ required by us is close to these values. 
Hence, for this range, ti is almost purely right-chiral {6^ ~ 90°). No mixing is required 
for the purpose and hence At values are seen to be within a couple of hundred GeVs. 

The problem is now to understand if the imprints of these branching fractions can be 
read out by studying suitable final states at the LHC. As indicated earlier, this is going 
to be rather challenging given that the third generation squarks all lead to very similar 
final states. Thus, disentangling individual contributions (read, contaminations), which is 
so necessary to unravel the sector, is a complicated proposition. This is more so since, in a 
bottom-rich environment, identifying multiple bottom quarks could very well hold the key. 

To simulate SUSY events we use the event generator Pythia v6.420 |58]. Pythia is 
interfaced with the framework SUSYHIT ^] that in turn uses the SLHA |66j protocol 
to integrate Suspect 2.31 |67j, the popular SUSY mass-spectrum generator and SDECAY 
|68j and HDECAY [69] which calculate the decay branching fractions of the SUSY par- 
ticles and the Higgs bosons respectively. To simulate the SM backgrounds discussed in 
section 3, partonic events are generated with Alpgen v2.13 |70j We have used CTEQ6L 
|71j parametrisation of the parton distribution function (PDF) via Pythia's interface to 
LHAPDF v5.7 |72j . The renormalisation/factorisation scale is set to ^/S, for both signal 
and background analyses. In case of background processes in Alpgen for which ^/S is not 
available as a choice for the renormalisation/factorisation scale, the default option for the 
same is used. To attempt a somewhat realistic treatment of the final state objects like jets, 
leptons, photons and the missing transverse energy we interfaced AcerDET vl.O |73) as the 
fast detector simulator. 

Unweighted events for the signal processes from Pythia and that for the SM background 
processes from Alpgen are then showered using Pythia with initial and final state radiations 
on. To avoid possible double-counting, multijet events from the Matrix Element (ME) 
calculation in Alpgen are matched with the jets from the Parton Shower (PS) using the 
MLM prescription [74J, available in Alpgen as the default. In the entire exercise, a jet is 
defined with a minimum clustered energy in the calorimeter of 25 GeV having a cone size 
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of AR = 0.4 within |??| < 2.5. 

For a jet to be triggered as a 6-jet, it is required to have px > 5 GeV within \r]\ < 2.5 
and an isolation from a neighbouring jet of AR < 0.2 is required. Further, an average 
tagging efficiency of 50% is used for the 6-jets. 

Electrons and muons are triggered in AcerDET only if they have pT > 5 GeV, \r]\ < 2.4 
and have a lepton-jet separation AR^j > 0.4. Further, to ensure the purity of the candidates, 
energy deposited within a cone of radius AR < 0.2 about the candidate lepton was required 
to be within 10 GeV. 

We incorporate a generic set of kinematic cuts in our analysis following ATLAS |75j . 
The requirements for the final state objects and the employed cuts are as follows: 

• Two tagged 6-jets (with a ^-tagging efficiency of 50%), both with px > 40 GeV. 

• At least four jets with the hardest one requiring pj, > 100 GeV and the other three 
with pip > 50 GeV. On top of this, for an inclusive jet final state all jets should have 
pi^^ > 40 GeV. 

• For inclusive 1-lepton and SSDL final states, isolated leptons with p^ > 20 GeV are 
required. For OSDL final state a, both leptons are required to have p^ > 10 GeV. 
Also, for the SSDL and OSDL final states two and only two leptons are required. 

• missing > 150 GeV, 

On top the above set of cuts, we use cuts on two more kinematic variables which help 
reduce the SM background in an effective way. The first one is the so-called transverse 
mass of the system comprised of the lepton(s) and missing energy and defined as 

= ^{E'^ + - {pi + - (pi^ + ^^)2. 

Traditionally used to reconstruct a leptonically decaying 1^-boson, a suitable cut on Mt 
thus could efficiently reduce the VF-boson background from the SM processes. In the left 
panel of Figure |3] we show the Mt distributions for the two benchmark scenarios of Tables 
[T] and [2] and that for the SM background combined over all the contributing processes. 
The mx profile for the SM background (the profile in light-green) nose-dives beyond 100 
GeV with a subdominant tail extending up to 500 GeV. For the signal, the one for the 
high-mass case (in yellow) and the low-mass one (in blue) extend to 550 GeV and 650 GeV, 
respectively. As we can see, a cut of Mt > 200 GeV would efficiently eliminate the SM 
background where W bosons are associated. 

It is also known that the variable called effective mass (M^jf) where 

Meff = ^p'T + Y.PT+^T 

jets £eptons 

can also be a powerful generic discriminator in search of new physics with highly massive 
exotic states. In the right panel of Figure |3] we illustrate the Me// distributions for the 
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Figure 3: The transverse mass (left) and the effective mass (right) spectra for the summed-up back- 
ground (in light green/light grey) and the signal. The profile at the bottom in yellow/ash is for the 
signal in the high-mass benchmark scenario of Table [!} The profile in blue /black is for the signal in the 
low-mass benchmark scenario of Table [21 



signal and the SM background as described for the Mt curves and following the same colour 
convention. Studying the trends of the distributions in this plot, we require a flat minimum 
of 800 GeV for Mg// for analysing both high and low mass benchmark scenarios. In the 
latter case, it seems to be a little overkill though. However, we stick to this particular value 
to bring uniformity in the overall analysis (the SM background thus remaining the same) 
which indeed may render this part of the study only conservative. 

In the subsequent part of this section we demonstrate how the event counts for different 
final states vary not only as a function of 9^ but also when 9^ is varied while m^^ and m^^ 
are held fixed within ±10 GeV of the corresponding benchmark values. In all these cases 
we incorporated the kinematic cuts and the the 6-tagging efficiency as mentioned above. 

The variations in the mixing angles are effected by simultaneously tweaking w-g^ , , 
m^^ and At- The change in rates in different multi-lepton and jet final states reflect the 
sensitivity of the rates to varying mixing angles, masses being held fixed. For this, we again 
take the spectra of Tables [T] and [2] as the reference ones with somewhat higher and lower 
squark masses, respectively. 

In Table [5] we present the different event rates for the high-mass scenario of Table [l} 
Note that the benchmark scenario here is chosen in a way such that 9^ ~ 67° which means 
the mixing angle is halfway between the maximal (45°) and what corresponds to a purely 
right-handed ti (90°). To be precise, 9^ = 67° refers to about a 15% admixture of in 
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ti. The mixing in the sbottom sector remains practically fixed at a rather low angle (~ 2°) 
thus rendering bi to be almost purely left-chiral. 





BR[6i ^ iiW-] 


26 + > 4j 
> li 


26 + > 4j + 
+ OSDL 


26 + > 4j + 
+ SSDL 


25.1° , 4.7° 


89.4 % 


66 


11 


4 


46.6° , 46.0° 

^LR ^ ILR 


78.8 % 


74 
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45.5 % 
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From Table 1 
67.0° , 1.9° 

XLR IL 
''1 > "1 


64.0 % 


65 


14 


4 


84.7° , 2.6° 

fR IL 


9.4 % 


58 
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3 


87.9° , 87.2° 


0.0 % 
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4 






Table 5: Variations of events rates (for the spectrum in Table [T]) in different leptonic final 
states containing two tagged 6-jets along with other light quark jets and missing transverse 
energy with different combinations of Oj: and O^. The kinematic cuts used are as discussed in 
the text. The event-rates correspond to an accumulated luminosity of 300 fb~^ and an average 
6-tagging efficiency of 50%. The last row presents the summed-up events rates from various SM 
backgrounds (see text for details) for the respective final states under the same set of cuts. 

In the first column of Table [5j he mixing angles in the stop and the sbottom sectors are 
presented along with the corresponding chiral contents of ti and bi. The suffix LR on ti and 
bi indicates that these mass eigenstates have some significant L and R contaminations while 
either L or i? as a suffix indicates that these are either almost purely left- or right-chiral in 
nature, respectively. It is important to note the mixing angles in the first column in this 
Table. These combinations of angles are again obtained by varying t^Q^^ , ^fei? ' ^tn^ ^* 
Af, while holding m^^ and m^^ within ±10 GeV of the respective benchmark values. Same 
procedure is taken for subsequent Tables in this section. Although our intention has been 
to demonstrate the effect of different combinations of mixing angles in the sbottom and 
the stop sectors spanning over purely left-chiral pairs to right-chiral ones passing through 
intermediate situations, we cannot really realise a situation when both bi and ti are almost 
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pure left-chiral states and still allowing for the necessary splitting between m^^ and m^^. 
The reason is quite clear and is as follows. Since both bL and tL are states determined 
by Qs^, the mass-eigenstates bi and ii dominated respectively by them are bound to be 
rather degenerate (~ ). The only splitting between them is due to the SU(2) D- 

term (oc cos 2/3) originating from the different electroweak (isospin (T3), electric charge) 
quantum numbers the left-chiral states carry. The splitting between m^^ and (i-e., the 
one that ensures m^^ > mf_^+W) that we are particularly interested in, in the present study, 
is thus not achievable unless we allow for some mixing in the top squark sector. This is 
what is reflected in the first row of Table [s] The angle 9^ ~ 24.7° indicates that ii has some 
right-chiral admixture but still dominated (since, % < 45°) by the left-chiral state. This 
observation would have ramifications beyond the present context. It implies that given a 
mass-splitting between bi and ti, combination of arbitrarily small 0^ and 9^ (i.e., arbitrarily 
large left chiral components in both bi and ti) are not feasible. This is true even when we 
minimally supersymmetrise the SM without subjecting it to further constraints. This is a 
kind of entanglement we like to point out. In other words, the phenomenon rules out the 
possibility of having a full-strength coupling in the on-shell decay 61 — )• tiW (or, for that 
matter, ti — )• biW, for a reverse hierarchy of masses between ti and 61). 

From Table [5] it can be seen that the variations in the event counts of different leptonic 
final states as the combination of mixing angles vary are anything but drastic. This is 
not unexpected. The reason behind this is that when the branching fraction of 61 — )■ 
iiW decreases, the same in other decay modes, e.g., 61 — tx^^ and 61 — 6x2; start 
increasing (see Figure [l|. The decay width for bi — t- iiW could be directly affected by 
double suppression from both 9^ and 9^ at the same time while its other modes of decay 
only see 0^ directly. Thus, there can be a possible sharing of branching fractions among 
these modes all of which contribute to the leptonic final states discussed here. Hence, we 
may reasonably expect a less drastic variation of the event counts with varying 9^^ and 9^. 
when these modes are open. 

Thus, an appropriate setup to improve the sensitivity to the mixing angles and hence to 
see the tell-tale signatures of the variations of mixing angles on the event rates of various 
multilepton final states is to have a situation when the decay modes bi — >• txi^ and 
61 — )• 6x2 kinematically closed. This we realise by making M2 larger than m^^ (M2 = 750 
GeV) so that the lighter chargino and the second lightest neutralino becomes heavier than 
61. We also set jU = 1 TeV (changing it from 750 GeV, as was for the benchmark scenario of 
Table [s]) so that the compositions of the chargino and the neutralino sectors remain more 
or less unchanged. 

The resulting variations are presented in Table [6j Branching fractions are now only 
shared between the decay modes 61 — t- iiW and bi — t- 6x1 • Events from these two de- 
cay modes could, to a good extent, may unambiguously tell us about the respective decay 
branching fractions which, in turn, would be indicative of the mixing angles involved. How- 
ever, bi — )• bxi would lead to leptonically quiet events with 6-jets and f^x which are hard 
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to identify over the huge SM background, particularly damaging one being of QCD origin 
with 'fake' 6-jets. Thus, one needs to concentrate on how the rates in different multilepton 
final states are varying to get an idea about the mixing angle(s) in the sbottom and stop 
sectors. The combinations of angles are kept close to the corresponding ones of Table [5j 
As we can see, the branching fractions to tiW mode, shown in column 2, displays a much 
clearer variation with the angle-combinations when compared to Table [5} With increasing 
right-chiral component in both bi and ti, the branching fraction 61 tiW gets diminished 
straight-away. Note that this decay mode of sbottom is the only source of leptons from 
61 down the cascade. The other one being bi — )• 6x1 1 this would only lead to 6-jets and 
missing energy in the final state. Thus, a more drastic (compared to Table [5] decrease in 
the number of events in different leptonic final states is expected as right chiral components 
increase. This is exactly what Table [6] indicates. 
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Table 6: The same variations as shown in Table |5] but with nonuniversal gaugino masses where 
for bi, only the modes 61 iiW and bxi are open. Except for M2 = 750 GeV and fj, = 1 TeV 
all other parameters are as there in Table [T] 

In both Tables [5] and [6j the rates for both OSDL and SSDL final states are found to 
be on the lower side compared to the inclusive 1-lepton case. This is quite expected since 
final states with more number of leptons involve further suppression due to added leptonic 
branching. The combined SM backgrounds for the respective final states are indicated 
in the last rows of both the Tables. As can be gleaned from these numbers, except for 
the cases where both bi and ti have significant right-chiral components, the rates in the 
inclusive 1-lepton final state has significance above 5a at an accumulated luminosity of 300 
fb^^. On the other hand, only the OSDL rate with both ii and bi having maximal possible 
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left-chiral contamination (i.e., if^, bf combinations in the first row) passes the 5o" marl|^ 
at 300 fb~^ The signal rates in the SSDL mode are indeed low throughout. However, since 
the SM background for this can be virtually eliminated, seeing a few events would suffice. 
Note that, at this stage, even a two fold increase in the accumulated luminosity is not going 
to change the overall situation drastically, at least qualitatively, except for the fact that we 
may have then a handful of very clean SSDL events. The bottom-line of Table[6]is that there 
could be a clear imprint of the product of cos 6*^ and cos 9^ in the absolute and mutually 
relative rates of these leptonic final states, though at a somewhat high integrated luminosity. 
Thus, the corresponding mass value of the sbottom squark (around 700 GeV, along with 
the gluino of 1200 GeV) lives on the edge of explorability for such a study. This prompts 
us to take up an exactly similar study but now with a low-lying spectrum benchmarked in 
Table [2] and its non-universal (in terms of gaugino mass relationship) variant. 

In Tables [7] and [8] we present studies which exactly emulate the proceedings of Tables 
[5] and [6| respectively but using the spectrum of Table [2j With mj^^ « 500 GeV and 
rrig ~ 650 GeV, we definitely expect a larger yield of events in any given final state. That 
the numbers presented in both Tables [7] and [8] are of the same order as for their high- mass 
counterparts of Tables [5] and [6j reflects the fact that the former Tables have numbers for 
a much lower integrated luminosity of only 30 fb^^ as compared to 300 fb~^ for the latter 
ones. SM backgrounds also get scaled down by this luminosity factor. Tables [7] and [8] clearly 
demonstrate how the sensitivities to mixing angles can be better studied with lower masses 
of the involved SUSY particles for all the final states discussed. 

The results presented in Tables [7] and [8] have direct correspondences to their high-mass 
counterparts presented in Tables [5] and [6j Note that, in Tables [6] and |8j even when both 
ti and bi are both almost purely right chiral, there are still significant number of leptonic 
(in particular, inclusive one-lepton events) events present, which is not quite expected had 
hi ^ tiW been the sole source of leptons in the final state. In fact, in our simulation not 
only 61-pairs but also gg and gbi pairs are included, as pointed out in the beginning of 
section 3. Thus, most of the leptons, under such a circumstance, are coming from top- 
squarks produced in the decay of the gluino along with the ones from the decay of SM top 
quark obtained under the SUSY cascade. This brings us to an important issue. 

Of particular interest is the cascade 61 — )• ii W~ — ?• t Xi W~ ~^ ^ Xi W~ — )■ 
b + ieptons and/or jets + ^t- Here, note that the decay of a single 61 gives rise to 
two VF-bosons (with opposite charges): one coming directly from the decay of bi while 
the other appearing in the decay of a top quark (originating from a top squark) further 
down the cascade. Thus, if both bi-s undergo such a cascade, at some stage there would be 
four (4) W bosons there. Possibility of identifying more of them could provide us with a 
remarkable handle to remove the above-mentioned 'spurious' leptons that do not originate in 



^This assumes a Gaussian estimation of the significance. However, for some of these low yields a Poissonian 
treatment would be more appropriate. The main issue here, however, is to take a simple note of the depleting 
counts in the dilepton final states for an integrated luminosity of 300 fb~^. 
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Table 7: Variations of events rates (for the spectrum in Table |2]) in different leptonic final 
states containing two tagged 6-jets along with other light quark jets and missing transverse 
energy with different combinations of 6i and O^. The kinematic cuts used are as discussed in 
the text. The event-rates correspond to an accumulated luminosity of 30 fb~^ and an average 
6-tagging efficiency of 50%. The last row presents the summed-up events rates from various SM 
backgrounds (see text for details) for the respective final states under the same set of cuts. 

the cascades of bi. Even the model-background originating from direct (in pairs) production 
of ti, which has not been considered in this work, can be removed if multiple W-s can be 
reconstructed in the final states discussed here. This would facilitate direct probe to the 
coupling bitiW thus reflecting on the mixings in the bottom and the top squark sectors. 

To the best of our knowledge dedicated study in this line is still lacking. One can take 
useful cue from some recent studies on multi-top final states [TH [771 EH] . As pointed out 
in these works, reconstructing multiple top quarks is an inherently difficult proposition, 
particularly in the complex environment of the LHC where the final state objects over- 
lap. However, identifying multiple VF-bosons is not expected to be more complicated than 
tracking down multiple top quarks. This is because, identifying multiple top quarks involves 
successive reconstructions of first, W bosons and subsequently, the individual top quarks 
they are coming from. Also, note that with our choice of soft SUSY parameters, ti always 
decays to txi thus making one of the W bosons in the cascade always coming from a top 
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Table 8: The same variations as shown in Table 5 but with nonuniversal gaugino masses where 
for 6i decay, only the modes 6i tiW and bxi are open. Except for M2 = 550 GeV and /i = 700 
GeV, values of all other parameters are same as in Table |2} 

quark. This may offer some degree of simplicity while undertaking a feasibility-study of 
exploring such a cascade. It can be foreseen that capability of identifying these W bosons 
would prove to be crucial in a generic scenario where 61 can decay to other channels {t Xi^ , 
6x2) for example) which ultimately give rise to identical final states. 



5 Summary and Outlook 

Physics analyses of recent LHC data hint essentially to super- TeV squarks from the first two 
generations and to a somewhat massive gluino. However, till now, the LHC experiments 
are not as sensitive to searches for the sbottom and the stop squarks. Consequently, much 
lighter sbottom and stop squarks are still allowed. This is exactly under such a situation 
when searches for them assume special significance. 

In this work, we aim for a rather conservative approach. We presume a scenario where 
only the lightest of the bottom and top squarks have significant cross sections at the 14 
TeV run of the LHC. Gluinos are taken to be intermediate in mass such that they could 
contribute to our final states only in a moderate way. On the other hand, squarks from 
the first two generations are considered to be heavier than a TeV following recent analyses. 
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These are, however, stih within the reach of LHC with reasonable cross sections and are able 
to contribute to the final states considered by us in different ways. However, in the spirit 
of the present work, to be conservative, we neglected those contributions. The compulsion 
is then to learn from the limited, nonetheless crucial, piece of information offered through 
the squarks from the third generation. 

It is pointed out that requiring some splitting between the the lighter sbottom and 
the stop eigenstates may set in some kind of an entanglement between the two sectors. 
This is because of the common left-chiral soft mass that enters the diagonal terms of the 
mass-squared matrices for both the sectors. One possible fall out of such entanglement is 
that the mixing angle in one of these sectors may constrain the corresponding one in the 
other sector. The degree of such entanglement may depend upon quite a few factors. The 
absolute masses of bi and ti and the splitting between these masses are two such important 
ones. 

In this work we demonstrate the phenomenology of the lighter sbottom squark. The 
reference decay mode considered is bi — )• tiW. This decay mode is naturally enhanced 
when both hi and ti have significant left-chiral content in them. Sensitivity of the resulting 
phenomenology to variations of the chiral contents of bi and ti would thus help probe their 
compositions. We point out the situations under which a clean study is possible. Under 
involved situations, we suggest that our ability to identify more number of VF-bosons in 
the final state would hold the key. We also stress that observations in multiple final states 
naturally would facilitate understanding. We demonstrate the issue with somewhat heavy 
sbottom and stop followed by the lighter ones at the 14 TeV LHC run. For the heavier 
spectrum we required an integrated luminosity in the ballpark of 300 fb~^ while for the 
lighter spectrum 30 fb~^ may be good enough. 

Maximally (even moderately) mixed top and bottom squarks is something outside the 
realm of much popular CMSSM/mSUGRA frameworks. Any hint of such mixings would 
firmly indicate departure from these scenarios. As exemplified in Tables [6] and [8| it may 
happen, that such a study may crucially bank on the nonuniversal masses in the gaugino 
sector thus exposing another crucial piece of information on the SUSY spectrum in the 
same go. 

A further entanglement is envisaged between the chargino/neutralino sector and the 
sector comprising of the third generation squarks. Here, the common agents responsible 
are fj, and tan /3. The role of fi in such entanglements is expected to be tempered by 
the value of tan /3 and hence could well be limited to scenarios with light stop and sbot- 
tom. Further, due to the opposite roles played by tan/3 in determining the mixings in the 
stop and sbottom sectors, sbottom sector may see somewhat larger correlation with the 
chargino/neutralino sector. It is important to note that since both mixings and masses in 
the latter sector are controlled (to different extents) by fi and tan /3, the entanglements can 
easily have both kinematic and dynamical implications. Moreover, with squarks from the 
third generations involved, they respond to both gaugino and higgsino components of the 
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charginos/neutralinos in characteristic ways under SUSY cascades. 

In this work, we have not considered a substantial, positive higher order correction (at 
NLO or NLO+NLL combined) to the SUSY cross sections. Consideration of this, could 
increase the cross section significantly (by ~ 35%) for our benchmark scenarios. However, 
to be conservative in our approach, and that we are unable to consistently treat the SM 
background on a similar footing, we postpone this to a future work. Further, we have 
not included possible contributions from 62 and ^2 either from direct production or from 
cascades. The former can very well be small while contribution from the cascades may still 
be appreciable under favourable circumstances. Neither do we discuss a possible reverse 
hierarchy of m^^ > m^^ where a similar phenomenology with stress on the decay ti ^ biW 
can be studied. Also, a dedicated analysis could have been undertaken for the 7-TeV run of 
the LHC. Apriori, given that the cross sections for the processes relevant for such a study 
are much smaller at 7 TeV, only some legitimate low-lying spectra might be of interest. 
However, even in that case, the requirement for total integrated luminosity would be in the 
ballpark of < 100 fb^^ which the LHC is not foreseeing for its 7- TeV run. 
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